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Arising from emergent electromagnetic field of magnetic skyrmions due to their nontrivial topology, the 
skyrmion Hall effect might be a roadblock for practical applications since any longitudinal motions of 
skyrmions in nanotrack is accompanied by a transverse motion. A direct consequence of such an effect 
is easy destruction of skyrmions at the nanotrack edges during their fast motions along the nanotrack, 
despite their topological protection. Here we propose an entirely novel solution of completely inhibit¬ 
ing such skyrmion Hall effect without affecting its topological properties based on a antiferromagnetic¬ 
coupling bilayer system. We show that a pair of magnetic skyrmions can be nucleated in such a bilayer 
system through vertical current injection or converted from a current-driven domain-wall pair. Once 
nucleated, the skyrmion pair can be displaced through current-induced spin torque either from a vertical 
injected current or in-plane current. The skyrmion Hall effect is completely suppressed due to the cancel¬ 
lation of back-action forces acting on each individual skyrmion, resulting in a straight and fast motion of 
skyrmions along the current direction. This proposal will be of fundamental interests by introducing the 
bilayer degree of freedom into the system. Moreover, it provides an easy way to engineer the transport 
properties of the skyrmionic devices to achieve desired performance, making it highly promising for prac¬ 
tical applications such as ultradense memory and information-processing devices based on skyrmions. 


Since the first experimental observations of magnetic 
skyrmion lattices in bulk non-centrosymmetric magnets^^and 
filmP, there has been tremendous interest in these topolog¬ 
ically protected spin configurations with a quantized topo¬ 
logical numbeP^. However, the creation and transmission 
of isolated magnetic skyrmion in magnetic thin films are re¬ 
quired for any practical applications such as encoding in¬ 
formation in individual skyrmion to allow entirely novel de¬ 
vices and circuitr}32]. Significant efforts and progress have 
been made towards realizing such ultrathin film based on per¬ 
pendicularly magnetized magnetic layer | heavy metal struc¬ 
ture to h ost a rich variety of chiral spin textures including 
skyrmion j^^ l ^^ i The strong spin orbit coupling at the inter¬ 
faces between the magnetic layer with perpendicular magnetic 
anisotropy (PMA) and the underlying heavy metal layer pro¬ 
vides a sizeable D zyal oshinsky-Moriya interaction (DMI) to 
stabilize skyrmionJ®®^. 

To move a skyrmion in nanotrack for information¬ 
processing applications, a convenient and efficient way is by 
means of spin current which can transfer the angular momen¬ 
tum from itinerant co nduct ion electrons to the magnetic mo¬ 
ments of the skyrmioil®2^ However, one major roadblock to 
the manipulation and transmission of skyrmions in nanotrack 
is the skyrmion Hall effect, i.e., skyrmions exhibit the Hall ef¬ 
fect driven by spin currents due to the presence of the Magnus 
force which in turn originates from its nontrivial topolog})^. 
Thus a skyrmion will not move parallel with the direction of 
the current. Instead it will gain a transverse velocity with a 
magnitude proportional to the spin current density which dis¬ 
places the skyrmion towards the edge of the nanotrack. There¬ 
fore a skyrmion will be easily destroyed for a distance of 
much less than Iptm for a nanotrack made of typical magnetic 
layer I heavy metal systenP^^^. 

Since both the skyrmion Hall effect and its topological pro¬ 
tection arise from the same mechanism, i.e., its nontrivial 
topological number of ±1, it seems to be impossible to in¬ 


hibit the skyrmion Hall effect without breaking its topolog¬ 
ical protection. In this work, we propose a novel solution 
of two perpendicularly magnetized sublayers strongly cou¬ 
pled via the antiferromagnetic (AFM) exchange interaction 
with a heavy metal layer beneath the bottom magnetic layei^^ 
(Fig[^, b, c). When a skyrmion is created in one of the 
sublayers, another skyrmion is simultaneously created in the 
other sublayer under certain conditions. We refer to such a 
pair of magnetic skyrmions as a magnetic bilayer-skyrmion 
(Figj^i, e). Moreover, we show a bilayer-skyrmion can be 
displaced over arbitrarily long distances driven by spin cur¬ 
rents without touching the edges due to the absence of the 
skyrmion Hall effect, which is distinct from a skyrmion in 
the monolayer thin-film structure. In addition to address the 
above-mentioned dilemma, bilayer-skyrmions have many dis¬ 
tinct characteristics from conventional skyrmions in mono- 
layer thin film, which allow for versatile and multifunctional 
ultradense and ultrafast information processing and logic ap¬ 
plications. 

These features are true for both the current-in-plane (CIP) 
and current-perpendicular-to-plane (CPP) cases. The CIP im¬ 
plementation has been extensively studied in the past ten years 
for displacing domain walls or switching magnetization. Re¬ 
cently, the spin Hall effect has been demonstrated to be a 
more efficient means of manipulating magnetization. The cur¬ 
rent for the CIP case is along the nanowire axial direction. 
By contrast, the spin current is perpendicular to the heavy 
metal/ferromagnetic interface for the CPP case. In our nomen¬ 
clature for CPP, the spin current is perpendicular-to-film plane 
whereas the charge current can be either in the film plane 
(such as the spin Hall effect scenario) or perpendicular to the 
film plane (such as the perpendicular MRAM where a per¬ 
pendicularly magnetized polarizer is incorporated). The CPP 
scheme has been proven to be much more effici ent to m ove 
the domain wall or skyrmion than the CIP metho(P^^I®3. 

In this work, we first describe the nucleation process of an 
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FIG. 1: Schematics of the bilayer nanodisk, nanotrack and the bilayer-skyrmion. a, The bilayer nanodisk for the creation of skyrmions, of 
which the diameter is 100 nm. The spin current (polarized along —z) is injected into the top layer in the central circle region with a diameter 
of 40 nm. b, The bilayer nanotrack (500 nm x 50 nm x 3 nm) for the motion of skyrmion driven by current perpendicular to the plane 
(CPP). The charge current flows through the heavy metal along the x-direction, which leads to the generation of spin current (polarized along 
-hy) perpendicularly injected to the bottom layer due to the spin Hall effect. The skyrmion in the bottom layer is driven by the spin current, 
while the skyrmion in the top layer moves remotely due to the interlayer exchange coupling, c, The bilayer nanotrack (500 nm x 50 nm x 
3 nm) for the motion of skyrmion driven by in-plane current (CIP). The electrons flow towards the right in both the top and bottom layers, 
z.e., the charge currents flow along the -x-direction. The skyrmions in both the top and bottom layers are driven by the spin current. In all 
models, the thickness of both the top ferromagnetic (FM) layer, the bottom FM layer and the insulating spacer are equal to 1 nm. The top 
layer and bottom layer is antiferromagnetically exchange-coupled, where the initial state of the top layer is almost spin-up (pointing along 
+z) and that of bottom layer is almost spin-down (pointing along -z). d. Illustration of the bilayer-skyrmion in a nanodisk, which is a set of 
antiferromagnetically exchange-coupled skyrmions. e, Sideview of the bilayer-skyrmion along the diameter of d. The color scale represents 
the out-of-plane component of the magnetization, which is used throughout this paper. 


isolated skyrmion within the top layer by utilizing the CPP in¬ 
jection into a disk composed of two AFM-coupled magnetic 
sublayers. Owing to the interlayer AFM exchange coupling, 
another skyrmion will automatically emerge in the bottom 
layer. In so doing we explore various properties of a bilayer- 
skyrmion. Furthermore, we move a bilayer-skyrmion in the 
nanotracks either by the CPP or CIP. We also study the gener¬ 
ation of a domain-wall (DW) pair and the conversion process 
from a DW pair to a skyrmion in the bilayer nanotracks. 

Bilayer System coupled with AFM interaction 

Hamiltonian. We investigate the bilayer system where the 
top and bottom ferromagnetic (FM) layers are coupled anti¬ 
ferromagnetically by the exchange interaction, as illustrated 
in Fig[^. The Hamiltonian for each layer reads 

Ht=- Aintra m[ • mj + D • (m[ x mj) 

+ + i^DDb (1) 

where r is the layer index (r =T, B), represents the lo¬ 
cal magnetic moment orientation normalized as |m[| = 1 at 


the site i, and (i, j) runs over all the nearest neighbor sites in 
each layer. The first term represents the FM exchange inter¬ 
action with the FM exchange stiffness Aintra- The second term 
represents the DMI with the DMI vector D. The third term 
represents the PM A with the anisotropic constant K. iTooi 
represents the dipole-dipole interaction. There exists an AFM 
coupling between the top and bottom layers, 

-fainter — Winter ^ ^ '^i ' '^i • (2) 

i 

The sign of Ainter is negative refiecting that the interlayer in¬ 
teraction is antiferromagnetic. We assume that the spins in 
the top layer are pointing upward. Then the spins in the bot¬ 
tom layer are pointing downward due to the interlayer AFM 
couplings. 

Topological Number. The classical field rri^{x) is intro¬ 
duced for the spin texture in the FM system by considering 
the zero limit of the the lattice constant, a ^ 0. The ground- 
state spin textures are = (0,0,1) and = (0,0,-1). 
A magnetic skyrmion is a spin texture which has a quan¬ 
tized topological number. Spins swirl continuously around 
the core and approach the ground-state value asymptotically. 
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FIG. 2: Creation of skyrmions in the bilayer nanodisk and time evolution of the skyrmion number, a, Injection of skyrmions in the bilayer 
nanodisk (D = 4 mJ m“^) with/without interlayer AFM exchange coupling. A 0.5-ns-long spin current (P = 0.4) is injected into the top layer 
(denoted by solid yellow circles), followed by a 1.5-ns-long relaxation. The interlayer exchange constant Amter is set as 0 or -2 pJ m“^, 
whereas the corresponding interface exchange constant a equals to 0 or -2 mJ (see Supplementary Movies 1-2). b, The time evolution 
of the skyrmion number of the top and bottom layers in the nucleation process of skyrmions corresponding to a. c, Total micromagnetic 
energy Ptotai (including the intralayer exchange, interlayer exchange, dipolar, anisotropy and DMI energy) for a bilayer-skyrmion and the 
AFM-coupled ground state as a function of the DMI constant D. Relaxed state of (d) the top layer and (e) the bottom layer after the injection 
of a 0.5-ns-long spin current for various current density j and interlayer exchange constant Amter- 


The skyrmion is characterized by the topological number Qr 
in each layer, 

Qr = —^ [ d?xrri^{x) • {dxrri^{x) x dyrri^{x)) . (3) 
Ttt J 

We obtain Q^- = ±1 for a skyrmion in a sufficiently large 
area. We also call Q^- the skyrmion number. Even if the 
skyrmion spin texture is deformed, its skyrmion number does 
not change, as far as the boundary condition is not modified. 
It can be neither destroyed nor separated into pieces, i.e., a 
skyrmion is topologically protected. 

The spins in the top and bottom layers are tightly bounded 
due to the AFM coupling. Accordingly, if one skyrmion is 
created in the top layer, a second skyrmion is also created 
in the bottom layer (Fig|^) simultaneously. The topologi¬ 
cal number of the bottom layer is opposite to that of the top 
layer since all the spins are inverted, Qb = —Qt- See Fig|^ 
how these topological numbers evolve after the creation of a 
skyrmion in the top layer. 

Excitation Energy. We may determine numerically the 
spin profile of each skyrmion and estimate the excitation en¬ 
ergy based on the Hamiltonian iTtotai = Hj -\-Hb-\- i^inter- We 


compared the energy of one bilayer-skyrmion state and the en¬ 
ergy of the homogeneous state in Fig|^ as a function of the 
DMI. Similar result is obtained for any value of the interlayer 
AFM coupling strength. This is because the spin directions 
between the top and bottom layers are always opposite. 

The energy of the bilayer-skyrmion becomes lower than 
that of the uniform ground state for D > 4 mJ m“^, which is 
the threshold value. Consequently, the bilayer-skyrmion ex¬ 
citation is energetically favorable when the DMI strength is 
larger than the threshold value. However, spontaneous gen¬ 
eration of skyrmions does not occur due to the topological 
protection. Nevertheless, the topological protection can be vi¬ 
olated in condensed matter physics because of the existence 
of the lattice structure and the boundary of the sample. We 
may create skyrmions by leveraging these properties. 

Skyrmion Hall effect. It is well understoocP^ that the 
center-of-mass motion of a skyrmion is determined by the Fa- 
grangian. 


L = Lb-U, 


( 4 ) 
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FIG. 3: Creation of a bilayer DW pair, its conversion into a bilayer-skyrmion, and their motions driven by vertical current in a bilayer 
nanotrack. The length, wide width and narrow width of the nanotrack (D = 3.5 mJ m“^) equal 400 nm, 100 nm, and 20 nm, respectively, a, 
A local vertical spin current (j = 50 x 10^^ A P = 0.4, polarized along —z) is perpendicularly applied to the top layer of the narrow side 
(85 nm <x < 115 nm) before t = 50 ps. When the top and bottom layers are decoupled (Amter = 0 pJ m“^), only one DW pair is generated in 
the top layer. However, when the top and bottom layers are coupled (Ainter = -6 pJ m“^), a bilayer DW pair is created at t = 50 ps. b, A global 
vertical spin current is perpendicularly applied to the bottom layer (towards +z, polarized along +y) when the bilayer DW pair is created at t = 
50 ps. The current density j in wide part equals 5 x 10^^ A m“^, which is proportional to that inside the narrow part with respect to the ratio 
of narrow proportion (200 nm x 20 nm) to wide proportion (200 nm x 100 nm). The 45-dgree notches are employed to reduce the required 
current density for the DW-skyrmion conversion. When the global driving current is turned on at t = 50 ps, the local DW injection current is 
turned off at the same time. See Supplementary Movie 3. 


where is the Berry phase term of the spin texture, 

Lb = ^{XY-XY) = ^G-{RxR), (5) 

and U is the potential. Here, G = (0,0, G) is the gyromag- 
netic coupling constant representing the Magnus force with 
G = AttQ, and R = (X, Y) is the center-of-mass coordinate 
of the skyrmion. The Euler-Lagrange equation yields 

R = X F, (6) 

where Cz = (0, 0,1), and F = —VU is the force acting on the 
skyrmion. Consequently, a moving skyrmion feels the Mag¬ 
nus force and bends toward the direction perpendicular to the 
force. It is called the skyrmion Hall effect. 

The direction of the Magnus force is opposite when the 
sign of the skyrmion number Q is opposite. For instance, if a 
skyrmion is driven by the current along the Px direction it will 
be bended toward the Yy {—y) direction in the top (bottom) 
layer. 

Landau-Lifshitz-Gilbert-Slonczewski equation for CPP 

We may apply a CPP spin-polarized current injection from 
magnetic tunnel junctio n (MTJ) or the spin Hall effect in 
heavy metal layei®^!^. We numerically solve the Landau- 
Lifshitz-Gilbert-Slonczewski (LEGS) equation, which gov¬ 
erns the dynamics of the magnetization rrii at the lattice site 
i. By suppressing the layer index, it reads 


drrii 

dt 


I I Ti-pff dyrii 

- \-f\mi X iT- -f arrii x 
Y I7I u{mi X px nrii) - I7I u\mi x p), 


(7) 


where Hf^ = — Qi^totai/^rrii is the effective magnetic field 
induced by the Hamiltonian i^totai = Y Y TTinter, 7 
is the Gilbert gyromagnetic ratio, a is the Gilbert-damping 
coefficient originating from spin relaxation, u is the Slon- 
czewski STT coefficient, u' is the out-of-plane STT coef¬ 
ficient, and p represents the electron polarization direction. 
Here, u = | ^ | ^ 7 ^ with /io the vacuum magnetic permit¬ 
tivity, d the film thickness, Mg the saturation magnetization, 
and j the current density. We take — ^ direction for creating 
the skyrmion, while Yy direction for moving the skyrmion. 
The STT is induced either by injection from a magn etic tun¬ 
nel junction polarizer or by the spin Hall effecP^^^. 

LLG equation for CIP n—i—1 

Alternatively we may apply a CIP injectiorP^^I^ to move 
skyrmions. We numerically solve the LLG equation, 


drrii 

dt 


= — |7|mi X Hf^ Y arrii x 


+ ^-inrioir) • V)mi 


drrii 

dt 


2eM, 

p|a^/3 

2eM, 


[rrii X {j{r) • V)mi], 


( 8 ) 


where (3 is the strength of the non-adiabatic torque and a is 
the lattice constant. 

Creation of a bilayer-skyrmion by vertical spin current. 

We employ a CPP injection with a circular geometry in a nan¬ 
odisk. The spin-polarized current (polarized along —2;) is in¬ 
jected into the top layer in the central circle region, as illus¬ 
trated in Fig[^. 

We demonstrate how the spin textures develop in Fig|^. 
The spins start to flip in both layers following the spin cur- 
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FIG. 4: Motion of skyrmions in the top and bottom layers of a bilayer nanotrack. Top-views of motion of skyrmions at selected interlayer 
exchange coupling and times driven by (a) the CPP injection and (b) the CIP injection (see Supplementary Movies 4-7). The parameters of 
the bilayer nanotrack are 500 nm x 50 nm x 3 nm, D = 3.5 mJ . The skyrmions are initially created by MTJ skyrmion injector placed 
on the top layer at x = 100 nm. For the CPP case, the spin current (P = 0.4) in the bottom layer is applied along +z but polarized along +y. 
The skyrmion in the bottom layer moves towards right driven by the spin current, while the skyrmion in the top layer remotely moves due to 
the interlayer coupling. For the CIP case, the skyrmions in both the top and bottom layers are driven by in-plane spin currents (P = 0.4). The 
velocities of skyrmions in the top and bottom layers as functions of current density j with (c) large interlayer exchange Winter = -6 pJ m“^ and 
(d) small interlayer exchange Winter = -0.06 pJ m“^. The cross symbol denotes the decoupling and destruction of skyrmions in the top and 
bottom layers due to large current density and small interlayer exchange coupling, where the velocities are calculated before the destruction of 
skyrmion. 


rent injection only in the top layer. When there is no inter¬ 
layer AFM coupling, a skyrmion is formed only in the top 
layer (see Supplementary Movie 1). By contrast, a skyrmion 
is formed also in the bottom layer upon the current injection 
in the presence of the interlayer AFM coupling (see Supple¬ 
mentary Movie 2). 

We show the evolution of the skyrmion number in Fig|^. 
It oscillates at the initial stage for t < 0.2 ns, and rapidly 
increases to 1. The skyrmion remains stable even when the 
current is switched off, demonstrating that it is topologically 
protected. During this process, the spins in the top and bot¬ 
tom layers are always anti-parallel. A comment is in order. 
The saturated skyrmion number is not exactly Q = 1 but 


Q = 0.89. This is due to the fact that there is a background 
skyrmion number which originates from the tilting the edge 
spins. It is Q = —0.09. Accordingly, the calibrated skyrmion 
number is Q = 0.98, which is almost unity. 

We present a nucleation phase diagram of a bilayer- 
skyrmion pair as a function of the current density and the in¬ 
terlayer AFM coupling in Fig|^ and Fig|^. When the mag¬ 
nitude of the injected current is strong enough, the bilayer- 
skyrmion is created. This is due to the fact that spin flip costs 
a certain energy. On the other hand, if the interlayer AFM 
coupling is too strong, the bilayer-skyrmion is suppressed due 
to the fact that the nucleation fleld and the coercivity increases 
with the interlayer AFM exchange, leading to a larger current 
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density for nucleation. 

Creation of a bilayer-skyrmion from a bilayer DW pair. 

A magnetic skyrmion can be created from a DW pair by us¬ 
ing a junction geometr}P^. In this scenario we first make a DW 
pair into a nanotrack of the top layer through the local CPP in¬ 
jection with — ^ direction. We show how the spins start to flip 
in the top layer and subsequently in the bottom layer driven by 
the AFM exchange force in Fig Then, the bilayer DW pair 
is shifted by applying CPP current, as shown in the process 
from t = 50 ps to f = 120 ps in Figj^. Here we consider 
the vertical injection of a spin current towards and polar¬ 
ized along -^y in the bottom layer. The CPP injection moves 
the bilayer DW in the rightward direction. When the bilayer 
DW arrives at the junction interface (t = 170 ps), both the end 
spins of the DW are pinned at the junction, whereas the central 
part of the DW continues to move due to STT in the wide part 
of the nanotrack. Therefore, the structure is deformed into a 
curved shape and a bilayer-skyrmion texture forms at t = 190 
ps (see Supplementary Movie 3). 

Current-driven motion of a bilayer-skyrmion in a nan¬ 
otrack. The magnetic bilayer-skyrmion can be displaced by 
the vertical spin-polarized current as in the case of the mag¬ 
netic skyrmion. We may employ the CPP injection or the CIP 
injection to drive a bilayer-skyrmion. In general, a moving 
skyrmion is easily destroyed by the sample edges due to the 
skyrmion Hall effect. Therefore, the maximum velocity of 
skyrmion in FM nanotrack is typically much less than 10^ m/s, 
limited by the edge confining forced of ^ j J. 

The skyrmion in the top layer follows the motion of the 
skyrmion in the bottom layer even when the current is not in¬ 
jected into the top layer. This is because that two skyrmions 
are bounded by the interlayer AFM coupling. There is no 
skyrmion Hall effect for a magnetic bilayer-skyrmion. This 
can be explained as follows. If there is no interlayer AFM 
coupling, a skyrmion in the top layer moves left-handed 
and the skyrmion in the bottom layer moves right-handed. 
However, when the AFM coupling is strong enough, two 
skyrmions are tightly bounded and the Magnus forces acting 
on the skyrmions between the top and bottom layers are ex¬ 
actly cancelled. Accordingly the bilayer-skyrmion will move 
straightly. This mechanism works both for the CPP and CIP 
cases (see Supplementary Movies 4-5). 

We show the relation between the magnitude of the injected 
current and the velocity in Fig|^ The velocity is proportional 
to the injected current density. For strong enough current, the 
bilayer-skyrmion is destroyed and split into two independent 
skyrmions. This is because that the skyrmion Hall effect in¬ 
creases as the current increases, which acts as the repulsive 
force between two skyrmions (see Supplementary Movies 6- 
7). 

With strong interlayer AFM exchange coupling, the cou¬ 
pled skyrmions move along the central line of the nanotrack 
at a high speed of a few hundred meters per second, without 
any transverse motion. However, with small interlayer AFM 
exchange coupling, the skyrmions in the top and bottom lay¬ 
ers will be decoupled due to the fast motion of the skyrmion in 
the bottom layer driven by large current. Once the skyrmions 
in the top and bottom layers are decoupled, the skyrmion Hall 


effect becomes active, leading to the destruction of skyrmions 
in the top and/or bottom layer by touching edge. At the same 
time, in the CPP case, the skyrmion in the top layer will stop 
motion. It can be seen that the coupled skyrmions driven by 
CIP current doesn’t require a fine tuning of damping and non- 
adiabatic torque coefficients. 

It is also worth noting that the interlayer AFM coupling 
does not produce a mass of the bilayer-skyrmion. When the 
driving current is suddenly turned off, the bilayer-skyrmion 
stops high-speed motion immediately (see Supplementary 
Movie 8). On the other hand, we also investigated the case 
where the DMI constant D is different between the top and 
bottom layers. It is found that the results of current-driven 
motion of the bilayer-skyrmion do not change much since the 
DMI only changes the radius of the skyrmion (see Supplemen¬ 
tary Movies 9-10). The massless property and its robustness 
make the bilayer-skyrmion an ideal candidate for practical ap¬ 
plications. 

Perspectives 

We have presented a novel solution of inhibiting the Hall 
effect of skyrmions without affecting their topological proper¬ 
ties, by exploring a new device made of antiferromagnetically 
exchange-coupled bilayer nanodisks and nanotracks. Com¬ 
pared with the mostly investigated skyrmion in the FM mono- 
layer system, the bilayer-skyrmion exhibits entirely distinct 
characteristics with regard to the current-transport behaviour 
and robustness. First, it can move strictly along the direc¬ 
tion of the spin current flow, which makes it more appeal¬ 
ing for motions in nanowires for ultradense memory appli¬ 
cations. This is in high contrast with the case of monolayer 
skyrmion, where the skyrmion information carrier can be eas¬ 
ily destroyed by the edges of the nanotracks. Second, it will 
be immune to magnetic field perturbations which might be 
generated externally or internally within the device circuitry 
since the net magnetic moment is zero. Third, by introducing 
the bilayer degree of freedom, the transport properties of the 
device can be engineered to achieve desired performance. For 
example, the skyrmion Hall effect can be intentionally sup¬ 
pressed or enhanced by tuning the magnetic properties of in¬ 
dividual layers. This newly proposed solution of transporting 
skyrmion information carrier for arbitrarily long distances at 
much enhanced velocity may be very appealing for versatile 
applications such as ultradense memory and information pro¬ 
cessing. Similar ideas can be extended to multilayer or su¬ 
perlattice where the skyrmions are strongly coupled to realize 
a better manipulation of skyrmions in nanotrack or extended 
thin films. 

Methods 

Modeling and simulation. The micromagnetic simula¬ 
tions are performed using the Object Oriented MicroMag- 
netic Framework (OOMMF) including the Dzyalosh inskii- 
Moriya interaction (DMI) extended modulJ^®^!^. The 
time-dependent magnetization dynamics is governed by 
the Landau-Lifshitz-Gilbert (LLG) equation including spin 
torqud^^^. The average energy density is a function of 
M, which contains the intralayer exchange, the interlayer ex¬ 
change, the anisotropy, the applied field (Zeeman), the demag- 
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netization and the DMI energy terms. For micromagnetic sim¬ 
ulations, the intrinsic magnetic parameters are adopted from 
RefsP® Gilbert damping coefficient a = 0.3 and the value 
for Gilbert gyromagnetic ratio is -2.21 lx 10^ m s“^. 
Saturation magnetization Ms = 580 kA m“^, intralayer ex¬ 
change stiffness A = 15 pJ m“^, DMI constant D = 0 ~ 6 mJ 
m“^ and perpendicular magnetic anisotropy (PMA) AT = 0.8 
MJ m“^ unless otherwise specified. The interlayer exchange 
coefficient Ainter is set from 0 to -10 pJ m“^, whereas the cor¬ 
responding interface exchange coefficient a equals from 0 to 
-10 mJ m“^ (a = Ainter / 1 nm), where denotes that the 
interface is antiferromagnetically coupled. The field-like out- 
of-plane STT coefficient u' is set to zero. All samples are 
discretized into cells of 2 nm x 2 nm x 1 nm in the simu¬ 
lation, which is sufficiently smaller than the typical exchange 


length 4.3 nm) and the skyrmion size to ensure the numer¬ 
ical accuracy. For all simulation of current-driven skyrmions 
reported throughout this paper, the skyrmions are firstly cre¬ 
ated at the designed spot of the nanotrack (x = 100 nm) by a 
local spin current perpendicular to the plane of the top layer. 
Then the system is relaxed to an energy minimum state with¬ 
out applying any current. Next, we start the timer and the spin 
current (P = 0.4) is injected into the nanotrack with the geom¬ 
etry of current-in-plane (CIP) or current-out-of-plane (CPP) 
as shown in Figj^ In the configuration of CIP, the electrons 
flow toward the right in both the top and bottom layers, i.e., 
the currents flow toward the left, while in the configuration 
of CPP, the electrons fiow toward the top only in the bottom 
layer. 
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